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Abstract 

We present an approximate analytic solution of the Klein-Gordon equation in the presence of 
equal scalar and vector generalized deformed hyperbolic potential functions by means of param- 
eteric generalization of the Nikiforov-Uvarov method. We obtain the approximate bound state 
rotational-vibrational (ro-vibrational) energy levels and the corresponding normalized wave func- 
tions expressed in terms of the Jacobi polynomial Pn (x), where [A > — l,u>— 1 and x S [—1, +1] 
for a spin-zero particle in a closed form. Special cases are studied including the non-relativistic 
solutions obtained by appropriate choice of parameters and also the s-wave solutions. 
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I. INTRODUCTION 



In nuclear, molecular and high energy physics [1,2], one of the interesting problems is to 
obtain exact solutions of the relativistic wave equations like Klein-Gordon (KG), Dirac and 
Salpeter wave equations for mixed vector and scalar potential. The KG equation has also 
been used to understand the motion of a spin-zero particle in large class of potentials using 
different methods. This allows us to introduce two types of potential coupling, which are 
the four- vector potential (V) and the space-time scalar potential (S). 

For the case S = ±V, the solution of these wave equations with physical potentials 
has been studied recently. Furthermore, their exact solutions are possible only for certain 
central potentials such as Morse potential [3], Hulthen potential [4], Woods-Saxon potential 
[5], Poschl- Teller potential [6], reflectionless-type potential [7], pseudoharmonic oscillator [8], 
ring-shaped harmonic oscillator [9], Vq tanh 2 (r/r ) potential [10], five-parameter exponential 
potential [11], Rosen- Morse potential [12], and generalized symmetrical double-well potential 
[13], etc by using different methods. It is remarkable that in most of these works, the scalar 
and vector potentials are almost taken to be equal (i.e., S — V) [2,14]. Recently, interest 
in the solutions of the KG equation for the case where S(r) = ±V(r) has surged. It 
presents bounded solutions in the relativistic view, although the KG equation reduces to the 
Schrodinger-like problem in the nonrelativistic limit. However, the reduced equation can 
show the relativistic properties of the system. For the most recent contributions, one may 
consult the papers in [2,6,9,10] and references therein. Therefore, the choice S(r) = V(r) 
(positive energy states) produces a nontrivial nonrelativistic limit with a potential function 
E(r) = 2V(r), and not V(r). It represents the concept of the exact spin symmetry that occurs 
in nuclei {i.e., when the difference potential A(r) = V(r) —S(r) = 0). In the negative energy 
states (corresponding to S(r) = —V(r)) the nonrelativistic limit is the trivial interaction 
free- mode. It represents the concept of the pseudospin symmetry (i.e., occurs when the sum 
potential E(r) = V(r) + S(r) = 0) [15]. The spin symmetry is relevant for mesons [16]. 
The pseudospin symmetry concept has been applied to many systems in nuclear physics 
and related areas [15-17] and used to explain features of deformed nuclei [18], the super- 
deformation [19] and to establish an effective nuclear shell-model scheme [17,20]. This, of 
course, does not diminish the importance of such problems. It only limits it's contribution 
(with the proper physical interpretation) to the relativistic regime [14]. However, in some 
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few other cases, it is considered the case where the scalar potential is greater than the vector 
potential (in order to guarantee the existence of KG bound states) (i.e., S > V) [21-24]. 
Many authors have considered a more general transformation between the unequal vector 
and scalar potentials given by V(r) = Vq + /3S(r) (or S(r) ^ ±V(r)),where Vq and (3 being 
arbitrary constants of certain proportions have to be chosen after solving the problem under 
consideration (cf. Ref. [25] and references therein). Nonetheless, such physical potentials are 
very few. The bound state solutions for the last case is being obtained for the exponential 
potential with the s-wave (/ = 0) KG equation when the scalar potential is greater than the 
vector potential [21]. 

The problems connected with the molecular structure provide interesting and instructive 
applications of quantum mechanics, since molecules are considerably more complex in struc- 
ture than atoms. Two distinct problems arise in connection with molecular structure. The 
first is to obtain the electronic wave functions and potential energy functions of the nuclear 
coordinates. This problem can be solved analytically only in the simplest cases. The second 
is to obtain the solution of the nuclear motion equation. In this regard, the construction of a 
suitable potential function of a diatomic molecule is very important. It has been found that 
the potential-energy function for the lowest electronic states of actual diatomic molecules 
can be expressed by means of the Morse potential [26]: 



which has three adjustable positive parameters a, D and r e . At r = r e , it has a minimum 
value at zero and approaches D exponentially for large r. If ^ is somewhat smaller than r e , it 
becomes large (but not infinie) as r — > 0. This potential is important in the field of molecular 
physics describing the interaction between two atoms [26,27]. Overmore, progress has been 
made in the field of diatomic molecules and extensive use of the potential functions have 
been introduced [28,29]. At present the Morse potential is still one of the potential functions 
used most in molecular physics and quantum chemistry [30]. However, it has few asymptotic 
inaccuracies in the regions of small and large r. To avoid these inaccuracies, many works 
have been carried out in that direction to improve Morse potential [31]. In 1986, Schioberg 
[32] suggested hyperbolical (empirical) potential functions of the form: 




(1) 




(2) 
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where D, a and a are three adjustable positive parameters with D = D e /(1 — a) 2 (D e is 
the spectroscopic dissociation energy) and a ^ a/5, 5^0 is a parameter. In contrast to 
the everywhere-regular Morse potential, V + (r) is highly singular at the origin with 1/r 2 and 
1/r singularities. The two potentials behave similarly near the extremum point r = r e . It 
has the minimum value at the point 

r = r e = — arctan/i (o") ±1 , (3) 
a 

and approaches D exponentially for large r. Unlike the Morse potential (1), the empirical 
potential (EP) function V + (r) approaches infinity at the point r = 0. In the region of large r, 
it is closer to the experimental Rydberg-Klein-Rees (RKR) curve than the Morse potential 
for some diatomic molecules. For a diatomic molecular model instead of the nuclear model, 
we consider the reduced mass definition. If the nuclei have masses m\ and m 2 , the reduced 
mass is defined as fj, — mim 2 / (mi + m 2 ) and at this point the diatomic molecular model can 
be included to the the pseudospin symmetry concept. In this context, it is worth noting that 
the energy of the nuclear motion (in MeV) is widely separated from the energy in atomic 
vibration and rotation (in KeV). That is why, one is able to separate the two motions and 
study the atomic vibration and rotation separately from the nuclear motion. 

Lu et al [33] obtained an approximate solution for the Schrodinger equation of diatomic 
molecule oscillator with the positive sign EP functions given in Eq. (2) (5 = 1 case) by 
means of the hypergeometric series method. Further, rigorous energy eigenvalues and eigen- 
functions for the ID Schrodinger eqauation are also obtained using a similar method. Since 
there are no exact analytical solutions for the EP functions (I ^ 0) without an approximation 
to the centrifugal term ~ 1/r 2 , some approximation [34] was used to obtain these solutions. 
One of these approximations have been employed in solving the rotating Morse potential for 
any /-state [34,35]. The ro- vibrating energy eigenvalues of the EP functions were determined 
with a semiclassical (SC) procedure (the Bohr-Sommerfeld quantization condition) and a 
quantum-mechanical (QM) method (the Schrodinger equation) [33]. Furthermore, several 
approximation schemes have been developed to find better analytical formulas for Eq. (2) 
with 5=1 [33,36]. Overmore, Jia et al [37] have used the basic concepts of the super- 
symmetric quantum mechanics formalism and the functional analysis method to investigate 
approximately the pseudospin symmetric solutions of the Dirac equation for the arbitrary 
pseudo-orbital angular momentum number I and to obtain the bound state solutions for the 
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nuclei in the relativistic EP as a diatomic molecular model. Very recently, the NU method 
[38-40] was applied to solve the radial Schrodinger wave equation with the EP functions 
for / 7^ case [41]. The analytic solution is used to obtain the ro- vibrating energy states 
for selected H 2 and Ar 2 diatomic molecules using the relevant potential parameters and 
spectroscpic constants given in Ref. [32]. 

Over the past years, the quantum deformation [42] has been the subject of interest because 
of its relevance with applications in nuclei [43-45], statistical quantum theory, string beam 
theory and conformal field theory [46-49]. Recently, some authors have introduced few 
potentials in terms of hyperbolic functions [50,51] in the view of g-deformation [52]. 

Encouraged with the high performance of the above inter-molecular potential, we write 
the EP functions as 

e a(r-r e ) _|_ e ~a{r-r e ) ]_ _|_ g-2a(r-r e ) j _|_ g e ~2ar 

[coth a(r - r e )\ = - = - - e _ 2a(r _ re) = 1 - ge _ 2ar = [coth g (ar)} , 

(4a) 

V±(r,q) = D {1 - a [cotha(r - r e )] ±l } 2 = D jl - a 

where q = e 2are , giving a magnitude for q that is larger than one. For an inverse transfor- 
mation q = e~ 2aTe , the magnitude for q varies between zero and one. In the context of the 
quantum deformation [52] , the above form is similar to a g-deformed (perturbed) generalized 
deformed empirical potential (g-DEP/GDEP) functions. It is worth noting that the range 
of parameter q was taken as q > in [50] and has been extended to— l<g<0org>0 
or even complex by Ref. [51]. Such g-deformed potential functions have been introduced 
for the first time by Arai [42] for real q values. When q is complex, these functions are 
called the generalized deformed potential functions. In this paper, we intend to find the 
analytic solution of the KG equation for the equal scalar and vector g-DEP/GDEP with 
any orbital angular quantum number (1^0). The specific choice of V(r) = S(r) allows one 
to make KG equation approximately soluble for it's relativistic energy eigenvalues and wave 
functions. Further, it opens up a new approach of generating the non-relativistic solution 
which is found to coincide with the previous Schrodinger solution of Eq. (4). In the present 
calculations, we apply a parameteric generalization procedures of the NU method which are 
making our calculations straightforward and simple. 

The present paper is organized as follows. In sect. 2, we present a parameteric gener- 
alization of the NU method which holds for the exponential-type potentials. In sect. 3, 



1 ± qe 
1 =F qe 



—2ar 



-2ar 



(4b) 
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we obtain an approximate analytic NU bound state solution of the (3 + l)-dimensional KG 
equation for equal scalar and vector g-DEP/GDEP functions with arbitrary /-states. In 
sect. 4, we discuss two special cases, the vibrational (/ = 0) and the non-relativistic limit 
(Schrodinger solution). In Sect. 5, we calculate the ro- vibrating energy states for selected 
H 2 and Ar 2 diatomic molecules in the non-relativistic approach. Section 6 contains the 
relevant conclusions. 



II. NU METHOD 

The NU method is briefly outlined here and the details can be found in [38] . This method 
was proposed to solve the second-order differential equation of the hypergeometric-type: 

a\z)g'\z) + a(z)r(z)g'(z) + a(z)g(z) = 0, (5) 

where a(z) and a(z) are at most second-degree polynomials and r(s) is a first-degree poly- 
nomial. The primes denote derivatives with respect to z. To find a particular solution of 
Eq. (5), one can decompose the wave functions, g n i(z) as follows: 

g(z) = (f)(z)y n (z), (6) 

leading to recast (5) in the hypergeometric-type equation 

a(z)y f >(z) + T(z)y' n (z) + \y n (z)=0, (7) 

where 

\ = k + n'(z), (8) 

and y n i(z) satisfies the Rodrigues relation 

A d n 

In the above equation, A n is a constant related to the normalization and p(z) is the weight 
function satisfying the condition 

[a{z)p{z)\ = r(z)p(z), (10) 

with 

T (z) =t(z) +2tt(z),t'(z) < 0. (11) 



The weight function should be carefully chosen because it has an influence on the per- 
formance of orthogonal wave functions, orthogonal in the interval [0,1], of the type La- 
guerre Ln\z) and Jacobi P%*'^\z) polynomials etc. It is defined with a compact sup- 
port often called domain of influence which can be spheres in three-dimensions. Gener- 
ally speaking, the weight function commonly used is exponential function. Furthermore, 
the weighted integral and weighted average are defined by h(z) = J f{z)p{z)dz and 

g(z) = J f(z)p(z)dz/ J p(z)dz, if f(z) : O e E, real-valued orthogonal polynomial func- 

n n 

tions, respectively, with f(z) = a(z) and g(z) = pjf'^\z), Ln\z), ■ ■ ■ , etc. Since p{z) > 
and a(z) > 0, the derivative of t{z) needs to be negative [38] which is the essential condi- 
tion in making the choice of particular solution relevant to the real bound state solution. 
The other part of the wave functions in Eq. (6) is mainly the solution of the logarithmic 
derivative: 

(f)'(z) n(z) 



(f)(z) a(zY 

where 



(12) 



n(z) = \ [a\z)-r{z)\ ± sj\ W'(z) -r(z)} 2 - a(z) + ka(z). (13) 

is a polynomial of order one. The determination of k is the essential point in the calculation 
of 7i(z), for which the discriminant of the square root in the last equation is set to zero. This 
gives the polynomial ir(z) which is dependent on the transformation function z(r). Also, the 
parameter A defined in Eq. (8) takes the form 

A = A n = -nr'(z) - -n(n- 1) a"(z), n = 0, 1, 2, • • • . (14) 

In this regard, we can construct a parameteric generalization of the NU method valid for any 
central and non-central exponential-type potentials. We begin by comparing the following 
general hypergeometric equation 

[z (1 - c 3 z)] 2 g"{z) + [z(l- c 3 z) (ci - c 2 z)\ g'(z) + (-B lZ 2 + B 2 z - B 3 ) g{z) = 0, (15) 

with it's counterpart equation (5) to obtain 

t(z) =ci- c 2 z, 

a(z) = z(l - c 3 z) , 



a(z) = -B lZ 2 + B 2 z - B 3 . (16) 

where the parameters q and Bi (i = 1,2,3) are constants to be determined during the 
solution process. Thus, following the method, we may obtain all the analytic polynomials 
and their relevant constants necessary for the solution of a radial wave equation. These 
analytic expressions are given in Appendix A. 



III. THE KG SOLUTION OF EQUAL SCALAR- VECTOR GDEP FUNCTIONS 

In relativistic quantum mechanics, we usually use the KG equation for describing a scalar 
particle, i.e., the spin-0 particle dynamics. The discussion of the relativistic behavior of spin- 
zero particles requires understanding the single particle spectrum and the exact solutions to 
the KG equation which are constructed by using the four- vector potential A a (A = 0, 1, 2, 3) 
and the scalar potential (S) . In order to simplify the analytic solution of the KG equation, 
the four-vector potential can be written as A A = (A , 0,0,0). The first component of the 
four- vector potential is represented by a vector potential (V), i.e., A = V. In this case, 
the motion of a relativistic spin-0 particle in a potential is described by the KG equation 
with the potentials V and S [1]. For S = V case [39], the (3 + 1) -dimensional KG equation 
is recasted to a Schrodinger-like equation and thereby the bound state solutions are easily 
obtained by using the well-known methods developed in nonrelativistic quantum mechanics 
[2]- 

Let us now consider the (3 + l)-dimensional time-independent KG equation describing a 
scalar particle (spin-0 particle) with Lorentz scalar S(r) and Lorentz vector V(r) potentials 
which takes the form [2,14,53] 

P 2 op - (V(r) - E R ) 2 + (S(r)+mc 2 ) 2 ] ^ KG {T) = 0, (17) 

where m and Er denote the reduced mass and relativistic binding energy of two interacting 
particles, respectively, with P op = —iHV is the momentum operator. It would be natural 
to scale the potential terms in Eq. (17) so that in the nonrelativistic limit the interaction 
potential becomes V(r), not 2V{r). We follow Alhaidari et al [14] to reduce the above 
equation to the form [54] 



V 2 + 



h 2 c 2 



i \ 2 f 1 x 2 



2 V(r)-E R \ - l-S(r)+mc 



*Pkg(?) = 0. (18) 



Thus, after making use of the equal scalar and vector GDEP functions (S±(r; q) = V±(r; q)), 
Eq. (18) recasts to 

1 



h 2 c 2 



„ 2 d 2 2 8 1 
v 2 = — + -— + — 



H H + V ± {r,q))] jil> KG (?) = 0, 
Md6\ Sm m) + sin 2 6 dip 2 



£4 

3=1 



(19) 



dr 2 r dr r z 

where a\ = mc 2 — E R , a\ = mc 2 + E R . It is woth noting that the solution of the (3 + 1)- 
dimensional KG equation can be reduced to the solution of the Schrodinger equation with 
the following appropriate choice of parameters: af — > —E^ R and a|//i 2 c 2 — > 2yu//i 2 , where 
// = m\mif (mi +m 2 ) is the reduced atomic mass for the diatomic molecular system [14,54]. 
In addition, we take the interaction potential as in (4) and decompose the total wave function 
*Pkg(~^)i with a given angular momentum /, as a product of a radial wave function Ri(r) = 
and the angular dependent spherical harmonic functions Yi m (r): [53-56] 



VW^) = 9 -^Y lm (r), 



(20) 



with angular momentum quantum numbers being I and m. This reduces Eq. (19) into the 
form 

n2 ,/, , n \t2„2" 



rf^r) 
dr 2 



1 



fr 2 c 2 



2 2 i 2 r-i 

a 1 Q; 2 + a 2 L> 



1(1 + i)h 2 c 2 



i - a ( : i — ^— <,(>■) o, <? / o, 

(21) 

where ^tr^ is the centrifugal potential and the boundary conditions g(0) = g(oo) = 
as we are dealing with bound-state solutions. Moreover, if I is not too large, the case of 
the vibrations of small amplitude about the minimum, we can then use the approximate 
expansion of the centrifugal potential near the minimum point r = r e as [33] 



~' ( ' +1) ^+^^ eip( : a r ) ,+A 3 



1 =F gexp(— 2qt) 



± exp(— 2ar) 
1 =F gexp(— 2ar) 



where 



A = l 

A 1 = ±2 [exp(2ar e ) =f 1] |3 
A 2 = [exp(2«r e ) T if 



1 =F exp(— 2ar e 



2ar P 



8ar e 

- — r - 3 - 2ar e 

1 =F exp(— 2ar e ) 



1 =F exp(— 2ar e 
2ar e 

1 =F exp(— 2ar e ) 
2ar P 



(3 + 2ar e ) 



3 + 2ar e - 



1 =F exp(— 2ar e 
2ar P 



4ar e 
1 =1= exp(— 2ar e ) 



(22) 

(23a) 
(23b) 
(23c) 



and higher-order terms are neglected. In fact, Eq. (22) is the approximate expansion of the 
centrifugal potential ^*t^ and is valid for all r r e , the minimum point of V±(r) since r 
is not singular there. However, the expansion is not valid near the singularity point r = 0. 
Overmore, it is a good approximation for small vibrations around the equilibrium separation 
r — r e . When I ^ 0, we have to use an approximation for the centrifugal term similar to the 
non-relativistic cases which is valid only for q — 1 value [33,53]. However, for s-waves, we 
remark that the problem can be solved exactly and the solution is valid for any deformation 
parameter q. At this point, it is important to mention that very similar expressions to the 
above expression for the energy states have also been found over the past years for the 
hyperbolical (exponential-type) potentials with 5 — * 1 in Eq. (2) for V + (r) (cf. Ikhdair and 
Sever in Ref. [36]). Very recently, a new improved approximation scheme [36,39] for the 
centrifugal potential term 1(1 + l)/r 2 was proposed which appears to be very different from 
the ones used by Refs. [33,53,54]. 

Putting z = ±exp(— 2ar) e (±1,0) for V±(r), and defining the parameters 

B 1 = q 2 K 2 nl + Sf - qQt - £ B 2 = 2qK 2 nl - Q h B 3 = K 2 nl , (24a) 

with 



K m = 76rAl a l D (! " ff ) 2 + Kl ^^^ o + a\al > 0, (25a) 




e 



+ (25b) 



a 2 h 2 c 2 4a 2 r 2 
Aq 2 a 2 Da 2 + 1<k1 + A 2 + q 2 a 2 h 2 c 2 > 0. (25c) 



e 



2ahc W r, 
we obtain the hypergeometric wave equation 

where g(z) = g(r). If we apply the previous transformations, the above expressions reduce 
into their non-relativistic limits: 



/v„, - K„, = ^\l2fiD (1 - a) 2 + l -^±^-A - 2fxE NR > 0, (27a) 




? e 



2fiqD 1(1 + 1) 

^( 1 -^) + 4ro-/^i> (27b) 



o; 2 /i 2 4a 2 r 2 
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Si - 5, = ^-J^Dai + l l±^L^l + ^2 > o (27c) 
2cm y r 2 g 2 

and also when the deformation parameter q = 1, the above equations reduce to their coun- 
terparts as in Refs. [33,41]. 

Now comparing Eq. (26) with Eq. (5), we obtain particular values for the set of constant 
parameters given in Section 2: 

ci = 1, c 2 = c 3 = q, c 4 = 0, c 5 = -|, c 6 = g 2 ^, + S 2 - 

c 7 = -2gK 2 , + Q h c$ = Kl h c 9 = Sf, 

2~ 

cio = 2if n , = 2c 12 > -1, en = -5/ = 2ci 3 - 1 > -1, 

Q 

Cia = > 0, Cia = ^ (<$ + §) > °- ( 28 ) 
Using Eqs. (28) together with Appendix A, we find the following particular physical solutions 
for the parameters: 

n(z) = K nl -(l + qK nl + S l )z, (29) 
k = -Qi - 2K nl S h (30) 

and 

t(z) = l + 2K nl -2(q + qK nl + £,) z, (31) 

where r'(z) = — _2 + qK nl + < which gives possible real solutions. In what 
follows, from Appendix A, we find the ro- vibrational energy equation with the aid of (28) as 

\ 2 

i 



1) - 


Qi 


1 






Q 


4 





+ n + 9 

2K nl = ^ » : ^ % J-, q ± 0, (32) 

f + n+\ 



which can be written more explicitly as 



2 X \ (mc 2 + E R ) D (1 - a) 2 + ^ + A + m 2 c 4 - £ 2 



4 (mc 2 + £*) Da + ^±ip - f ) - (e + aftc (2n + l))' 

E + ahc(2n+ 1) 



(33a) 



E = \ {mc 2 + Da 2 + l ( l + 1 ) H2c2 ^l + a 2^2 c 2 ? ( 33b ) 
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where n — 0, 1, 2, • • • and I — 0, 1, 2, • • • signify the usual vibrational and rotational angular 
momentum quantum numbers, respectively. 

Let us now turn to the calculations of the corresponding wave functions for the potential 
under consideration. Thus, referring to the general model in Appendix A, the explicit form 
of the weight function reads 

p{z) = z 2R - l (l-qz)^\ (34) 
which gives the first part of the wave functions (6) as 

y n (z) - pS Knl ' lqSl \l =F 2qz), K nl > 0, S t > 0, (35) 

with the essential requirement that 2K n i > —1 and -Si > — 1. For example, if q > then 
Si > — | and if q < then < Si < — |. Also, the second part can be found as 

1 

and, hence, the unnormalized wave functions are being expressed in terms of the Jacobi 
polynomials as 



) ^ z K nl{l T qrffr+i), Knl > 0, - (5, + |) > 0, (36) 



^(z) =Ar„^(l T g^)^ i+ ^pi 2i "" i '^ i) (l-2gz), z e[0,l/g] (37) 

where M n i being the normalization constants and P n " ,q (1 — 2qz) = 
^±^ 2 Fx(-w, 2K nZ + |Si + n + 1, 2^ ni + 1; gz) with (m) n = is Pochhammer's 

symbol. For example, if q > 1 then z£ [0, 1/q] and if g < — 1 then z G [l/g,0] lie within or 
on the boundary of the interval [—1, +1] . 

Hence, the total wave function of the g-DEP/GDEP functions is 

i)±(~r) =Af n i- [±exp(-2ar)] Xni [1 - ±q exp(-2ar)]^ i+ ^ 
r 

x Pn K nl ' lqSl \l - ±2qexp(-2ar))Y lm (r). (38) 
where the normalization constants Af n i are calculated explicitly in Appendix B. 



IV. DISCUSSIONS 

In this section, we are going to study two special cases of the energy eigenvalues given 
by Eq. (34). First, we consider the s-wave (I = 0) vibrational energy equation: 



2\j (mc 2 + E R ) D (1 - a) 2 + m 2 c 4 - E 
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4 (mc 2 + E R ) Da - (y4 (mc 2 + E R ) Da 2 + a 2 h 2 c 2 + ahc (2n + 1))' 



v/4 (mc 2 + Z)(; 2 + a 2 ft 2 c 2 + aftc (2n + 1) 



(39) 



where n — 0, 1, 2, 



: , where n max is the number of bound states for the whole bound 



spectrum near the continuous zone. n max is the largest integer which is less than or equal 
to the value of n that makes the right side of Eq. (39) vanish, that is, 



n 



ax 



-1 - 



AD 



a 2 h 2 c 2 



mc 2 + E R )a 2 + l + 



AD 

a 2 h 2 c 2 



(mc 2 + E R ) a 



(40a) 



(40b) 



If a = 1, then E, 



(R) 



mc 2 . The corresponding normalized wave functions can easily be 



found directly from Eq. (38) as 



R^\r) = N n - [±exp(-2ar)] fe " [1 - ±q exp(-2ar)]^ 



x Pn kn ' lqT> (l - ±2gexp(-2«r)), 



(41) 



where M n are the normalization constants and calculated in Appendix B. 

Second, we discuss the non-relativistic limit of the energy eigenvalues and wave functions 
in the non-relativistic limit. Obviously, the currently calculated KG solutions, under the 
previously mentioned transformations, can be reduced to their associated Schrodinger ones 
for the GDEP functions as 

(NR) _ n ,1(1 + l)ft 2 



E 



(R) 
nl 



E 



nl 



= D e + 



2/ir 2 



a 2 h 2 

2/i 

or more explicitly as 



(f) 2 -f-Hf+-+0 ; 



f + n+l 



-A) 



, n,/ = 0,l,2, 



(42) 



E 



(NR) 
nl 



D e+ l -^A 



2/ir, 



a 2 h 2 

2/i 



a 2 ° ~T~ ia 2 r 2 \ q 2 q J \"' ^ 2 ^ \ h 2 a 2 ° ^ ia 2 r 2 q 2 



2/iD 

h 2 



21 2 



n 



4a 2 r 2 q 4 



(43) 



which is identical to Eq. (28) of Ref. [33] if one sets 5 = 1 and q — 1. This represents 
the approximate Schrodinger solution of Eq. (2) for the ro-vibratinal molecules. The non- 
relativistic limits for the vibrational energy states (I = 0) read 
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and the corresponding unnormalized wave functions from Eq. (38) are 

V>±( "r ) = Mni^ [±exp(-2ar)]* nI [1 - ±q exp(-2ar)]^ i+ ^ 

X P { n K - ±2gexp(-2«r))yi m (f). (45) 

where and are defined in Eq. (27) and the condition for n max turns to become 

*£2> - A, (46) 

Thus, the finiteness of n max is reflected in the above condition. If a — 1, then n max — > 0. 

V. APPLICATIONS TO DIATOMIC MOLECULES 

We have calculated the non-relativistic energy states for the two selected H 2 and Ar 2 
diatomic molecules using energy equation (35) with q — ■> 1 and (23). The spectroscopic 
constants of these two molecules are given in Table 1. The vibrating ground state energy 
eigenvalues (in cm -1 ) for the H 2 molecule in the non-deformed EP functions V + (r) are 
found using the NU method for the potential parameters given in Table 2. Our numerical 
results obtained in the present NU model are listed together with the numerical results 
obtained by using SC (as Semi-Classical) procedure and a QM (as Quantum-Mechanical) 
method mentioned in Ref. [32] for various potential parameters. Obviously, as shown in 
Table 2, the results obtained in the present model are in high agreement with those obtained 
by QM. However, the SC procedure is proportionally different. Therefore, the differences 
between our results and SC procedure are less than 0.01 cm -1 , i.e., they are negligible 
because of these approximations: 1 a.m.u = 931.502 MeV/c 2 , 1 cm -1 = 1.23985 x 10~ 4 
eV and he = 1973.29 eV.A° [57]. The second application is applied to Ar 2 molecule. We 
confine our study to calculate the ro- vibrating energy states for the V+(r) potential using the 
following potential parameters: a = 25.23, 5 = 41.75 and a = 0.6604 (A )^ 1 [32] together 
with the parameters given in Table 1. For the previously given set of physical parameters, 

14 



n 



n r 



1 

2 



h 2 a 2 



a 2 + l + 



8fiD 
h 2 a 2 



a 



we plot the non-relativistic energy spectrum curve as a function of vibrational quantum 
number n as seen in Figure 1. Obviously, in a reference to Figure 1, the energy spectrum 
of the diatomic molecule Ar 2 approaches the value of D e as n approaches n max = 6.689 or 
n mffl = 6. This is also verified analytically from Eqs. (44) and (46). Moreover, a plot of 
the non-relativistic energy spectrum curve as a function of the potential strength D e for the 
above given set of physical parameters and n = for Ar 2 molecule is shown in Figure 2. 
The relationship is noticed to be nearly linear for any arbitrary value of vibrational quantum 
number n. The attractive energy value increases with the increasing potential strength. The 
splittings of the energy states of s-waves E + = E + (n ^ 0) — E + (n = 0) obtained by the NU 
method and SC procedures are presented in Table 3. The present results AE + (NU) from 
NU method and AE + (SC) obtained from the SC procedures are also compared with four- 
different experimental results labeled AE(a), AE(b), AE(c) and AE(d) taken from Ref. [32]. 
It is obvious from Table 3 that our results are very close with the experimentally determined 
values as well as the SC procedure results. Finally, the approximated rotating and vibrating 
energy states of the V + (r) given in Eq. (2) for the Ar 2 and H 2 molecules are also calculated 
for the I 7^ case. Table 4 shows the energy levels for vibrational (n = 0, 1,2,3,4,5) and 
rotational (I = 0, 1, 2) quantum numbers. 



VI. CONCLUSIONS 



To summarize, we have presented the approximate bound state energy eigenvalues and 
their corresponding normalized wave functions of the relativistic spin-0 particle in the radial 
(3 + l)-dimensional KG equation with equal scalar and vector g-DEP/GDEP functions by 
means of the parametric generalization of the NU method. We point out that the KG wave 
functions are found in terms of the Jacobi polynomials. The analytic expressions for the 
relativistic energy expression and the corresponding wave functions of this molecular system 
can be reduced to the well-known non-relativistic solutions and to the s-waves solutions as 
well. The relativistic energy E R defined implicitly by Eq. (33) is rather a transcendental 
equation and it has many solutions for any arbitrarily chosen values of usual quantum 
numbers n and I. The method presented in this paper is general and worth extending to 
the solution of other molecular interaction problems. The method is simple and useful in 
solving other complicated systems analytically without giving any restriction on the solution 
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of some quantum systems as is the case in the other models. We have also seen that for the 
nonrelativistic model, the approximate energy spectrum can be obtained either by directly 
solving the Schrodinger equation [41] or rather by even applying appropriate transformations 
to the relativistic solution as currently shown. We should emphasize that the approximate 
bound state energy spectrum obtained in the present work might have some interesting 
applications in different branches like atomic and molecular physics and quantum chemistry. 
The present solution is describing the inter-molecular structures and interactions in diatomic 
molecules [32-36,41,58,59]. The present study is also useful in calculating the vibrating 
energy for different radial n quantum numbers as well as the rotating energy for different 
orbital I quantum numbers. To conclude, the proposed g-deformation potential with a 
flexible and fixed value q (real or complex) can generate various potential models with 
various energy solutions. 
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APPENDIX A: PARAMETERIZED VERSION OF THE NU METHOD 

We complement the theoretical formulation of the NU method in presenting the essential 
polynomials, energy equation and wave functions together with their relevant constants as 
follows. 

(i) The key polynomials: 




(Al) 



k = - (c 7 + 2c 3 c 8 ) - 



(A2) 




(c 2 - 2c 5 ) z-2 [(y/c^ + c 3v /c^) z - y/c^\ , 



(A3) 




(A4) 
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(ii) The energy equation: 

(c 2 - c 3 ) n + c 3 n 2 - (2n + 1) c 5 + (2ra + 1) (v^ + c 3v ^) + c 7 + 2c 3 c 8 + 2^/^ = 0. (A5) 

(iii) The wave functions: 

p(z) = z cw (l-c 3 z) c '\ (A6) 

<j>(z) = z c ^(l - c 3 z) c ^, ci 2 > 0, cia > 0, (A7) 

j/ n (z) = Pi Cl0 ' Cll) (l - 2c 3 z), cio > -1, c n >-l,ze [0, l/c 3 ] , (A8) 

u(*) = AT„2 Cl2 (l - c 3 ^) Cl3 Pi Cll " Cll) (l - 2c 3 z), (A9) 

where the Jacobi polynomial P^ l ' u \x) is defined only for /j, > —1, v > — 1, and for the 
argument x G [—1, +1] and A/" n is a normalizing factor.. It can be expressed in terms of the 
hypergeometric function as 

PM(1 - 2s) = ^-i^ 2Fl (-n, 1 + ix + v + n; + 1; s) , (A10) 

where s G [0, 1] which lie within or on the boundary of the interval [—1, 1] . Also, the above 
wavefunctions can be expressed in terms of the hypergeometric function as 

u(z) = M n z Cl2 (l - c 3 z) Ci:i 2 Fi (-n, 1 + cio + cn + ra; c w + 1; c 3 z) , (All) 

where C12 > 0, Ci 3 > and z G [0, l/c 3 ] . 
(iv) The relevant constants: 

1 i 

c 4 = - (1 - ci) , c 5 = - (c 2 - 2c 3 ) , c 6 = C5 + Si, 

c 7 = 2c 4 c 5 - B 2 , c 8 = C4 + S 3 , c 9 = c 3 (c 7 + c 3 c 8 ) + c 6 , 

2 

cio = ci + 2c 4 + 2^ - 1 > -1, cn = 1 - ci - 2c 4 H ^ > -1, 

c 3 

C12 = c 4 + > 0, ci 3 = -c 4 + — (v^ - c 5 ) > 0. (A12) 

c 3 
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APPENDIX B: NORMALIZATION OF THE RADIAL WAVE FUNCTION 



In order to find the normalization factor Af n i, we start by writting the normalization 
condition: 

f 1 z 2K nl -i {l _ z ys l+ i \ P (^ nl ,2S l )^_ 2z )] 2 dz = i : (Bl) 
2a Jo L J 

where q — 1. Unfortunately, there is no formula available to calculate this key integration. 
Neveretheless, we can find the explicit normalization constant M n i- For this purpose, it is 
not difficult to obtain the results of the above integral by using the following formulas [59] 
i 



J (1 - s) M_1 s"- 1 2 F X (a, (3;r, as) dz 



r»r» 



? F 2 {v, a, (3- fx + v\ 7; a) , (B2) 



rox + i/) 3J 



oo 

and 2 F 1 (a, b; c; z) = V (a)T(b) ^T^ffi^ ff- Following Ref. [59], we calculate the normal- 



p=0 



ization constants: 



Til 



where 



T(2K nl + l)r(2S, + 2) ~ (-l) m (l + n + 2(K nl + 5,)) r(n + m) 

2 , 7 7 ~ ~ 



2«r(n) 



m=0 m! m 



+ 2X ni )ir (m + 2(^ + ^ + 1)) 



-1/2 



(B3) 



/„,= 3 F 2 (2X ni + m,-n,n + l + 2(X nJ + ^);m + 2(X ni + ^ + l);l + 2X ni ;l). (B4) 



Furthermore, the normalization constants for the s-wave can be also found as 

n -1/2 



r(2^„ + i)r(2?+ 2) ^ ( x + n + 2 (*» + ^)j( n + m ) 



2arfn) 



o m! m + 2/c n !r (m + 2(k n + s + 1) 



0r. 



, (B5) 



where 



= 3^2 ( 2k n + m, -n, n + 1 + 2(/c n + 3); m + 2(A; n + s + 1); 1 + 2/c n ; l) , (B6) 



and 



k n 



2ah 



mc 2 + E R ) D(l- a) 2 + m 2 c A - E 2 W 



2ah 



y/A (mc 2 + E R ) Da 2 + a 2 h 2 , n = 0, 1, 2, 



(B7) 



where E R is the solution of the transcendental equation (33). 
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FIG. 1: A plot of the non-relativistic energy spectrum curve as a function of the vibrational 
quantum number n for a given set of physical parameters for Ar^ molecule. 



FIG. 2: A plot of the non-relativistic energy spectrum curve as a function of the potential strength 
D e for a given set of physical parameters and vibrational ground state n = for Ar2 molecule. 



TABLE I: The spectroscopic constants of the EP for H2 and Ar2 molecules [26]. 

Parameters H2 Ar2 
D e (cm- 1 ) 38281 99.55 
r e (A°) 0.7414 3.759 
H (a.m.u) 0.50407 19.9812 
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TABLE II: The EP parameters of the V+(r) and the ground state energy, E^ (in cm 1 ) of the H2 
molecule. 

a S a {A )- 1 E + {SC) E + (QM) Present 

426.826 463.102 0.9327 2167.68 2168.93 2168.68 

47.294 102.341 0.6146 2153.69 2164.83 2164.45 

28.685 117.121 0.3826 2139.57 2157.69 2157.53 

21.250 213.212 0.1762 2124.29 2148.40 2147.53 



TABLE III: Comparisons of experimentally calculated s-states energy transition values 
A£'„ > o(cm _1 ) for n / — > n = together with the results of the SC procedure and the present 
NU method for the Ar2 molecule. 



n Present AE(a) AE(b) AE(c) AE(d) AE + (SC) 



1 25.808 


25.74 


25.49 


25.21 


25.56 


25.75 


2 46.079 


46.15 


45.63 


45.02 


46.00 


46.01 


3 61.472 


61.75 


60.70 


60.04 


61.32 


61.42 


4 72.536 


72.66 


71.33 


70.92 


71.52 


72.52 


5 79.733 


79.44 








79.79 


6 83.453 










83.59 


7 84.026 
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TABLE IV: Energy levels E n i(cm 1 ) for Ar 2 and H 2 molecules in V+(r) using the NU method. 
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